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q«^~ . The spectral properties of the Laplacian operator on "small- world" lattices, that is mixtures of 

' unidimensional chains and random graphs structures are investigated numerically and analytically. 

, A transfer matrix formalism including a self-consistent potential a la Edwards is introduced. In the 

t-H ■ extended region of the spectrum, an effective medium calculation provides the density of states and 

» | ' pseudo relations of dispersion for the eigenmodes in close agreement with the simulations. Local- 

ization effects, which are due to connectivity fluctuations of the sites are shown to be quantitatively 
^5 ■ described by the single defect approximation recently introduced for random graphs. 
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Ch ' I. INTRODUCTION 

P3 : 

Diffusion in random media is an an issue of great practical and theoretical importance which has a received a lot of 
attention in the past years jl]-^). The presence of disorder, e.g. impurities or random site energies, lattice defects, ... 
| may induce the localization of eigenstates and therefore have dramatic consequences on transport properties (|||||,(7| . 
The understanding of the localization transition, based on scaling theory for finite-dimensional systems has also 
greatly benefited from the study of models on tree- like architecture Jlfj , that is without any underlying geometry. 
Diffusion on such amorphous structures has also attracted attention as a toy-model of dynamical evolution in the 
phase space of complex, e.g. glassy systems [ fil] , fl2| . 

Recently, it was pointed out that mixed lattices, built from a mixture of finite-dimensional structures and random 
£j ' graphs could be of interest to modelize various real situations taking place in different scientific fields ]l3| . As a main 
attractive feature, these "small-world" models simultaneously combine short- and long-range interactions and could 
display some new and interesting properties stemming from this unusual coexistence [|13|-^5[ . 

From a physical point of view, this situation is already encountered in the field of polymers or more generally long 
■ molecules |16[] . Couplings between successive monomers give rise to a unidimensional system while (e.g. hydrody- 
namics) interactions between constituents in the three-dimensional space are of long range nature along the polymeric 
chain. Of course, the long-range links formed this way are strongly correlated and cannot be considered to be inde- 
pendent as in the small-world paradigm. Some successfull approximation schemes developed in the context of polymer 
theory, e.g. the self-consistent potential of Edwards |17 R| method precisely amount to neglect such correlations and 



nevertheless obtain reliable results JT(|. Following Zimm's wor k fl9[|, many studies have concentrated on the influence 



of such interactions on the relaxation dynamics of polymers |16| , |20| . However, these calculations can in a sense be 
considered as annealed since the hydrodynamical interactions tensor is averaged over chain positions before normal 
modes are computed and not later as it should. Such a pre-averaging approximation |2fJ may lead to qualitatively 
erroneous predictions especially for large eigenvalues and must be considered with caution. Diffusion on the small- 
world lattice can therefore be viewed as an elementary step toward the understanding of small times vibrations or 
equivalently the instantaneous normal modes pl|-p3[ of a polymeric chain around a quenched configuration. 

From a theoretical point of view, the small-world structure is also worth being investigated since it mixes two 
structures which can be studied in great mathematical details [ fl3| |. Diffusion on a one-dimensional lattice obsviously 
reduces to the study of dispersion relations for plane waves ||. As for the random graph structure, a large attention 
has also been paid to its spectral properties p^ , p^ -p8[. The main physical feature is the emergence of localized 
states centered on the geometrical defects of the graph, that is sites with abnormally low or large connectivities with 
respect to the average coordination number [ p6[ . As we shall see, the combination of both structures can be precisely 
studied and allows to answer some interesting questions concerning localization and the existence of pseudo-dispersion 
relations for the eigenmodes in presence of disorder ||^9| . 

This paper is organized as follows. In Section II, the small-world model is defined and related to its building bricks, 
that is a random graph superposed to a one-dimensional lattice. The main features of diffusion on the latters are briefly 
recalled and used to guess some spectral properties of the small-world Laplacian. Section III exposes the numerical 
results obtained from exact diagonalizations. The density of states as well as some characteritics of the eigenmodes, e.g. 
localization properties and autocorrelation functions are analyzed. We present in Section IV the analytical approach 
to attack the problem, based on a self-consistent transfer matrix formalism [ fiT|]l§f| . The region of the spectrum 
corresponding to extended eigenstates is studied in Section V by means of an effective medium approximation [fll,E9| . 
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Localization properties are unravelled in Section VI using a more refined scheme of approximation exploiting the 
geometrical defect mechanism exposed above Et| . 



II. THE SMALL- WORLD LATTICE AND ITS STRUCTURE 



A. Presentation of the model 



We consider N points Ai, i = 1, . . . , N on a one-dimensional ring 1Z. Each point Ai is connected to its 2K nearest 
neighbors Ai±i, Ai±2, . . . , Ai±K (where site labels are defined modulo N to ensure periodicity on 7Z). In addition, we 
superpose to 1Z a random graph Q with mean connectivity q defined as follows. Each pair of points (Ai,Aj) (such 
that the distance \i — j\ on the ring is larger than K) is connected with probability q/N and left unchanged with 
probability 1 — q/N. 

The small-world lattice S is the union of both graphs 1Z and Q. In average, S includes (K + q/2)N edges to largest 
order in N. The coordination degree C; of Ai is a random variable bounded from below by 2K. More precisely, Ci — 2K 
obeys a Poisson law of parameter q. This "small-world" lattice differs from the definition of Jll| but maintains the 
coexistence between short- and long-range links. In addition it leads to simpler analytical calculations than the model 
exposed in fl3|| . 

We then consider the Laplacian operator Wfj on S. For i ^ j, W$ = —1/2K if Ai and Aj are connected, 
otherwise. Diagonal elements, W§ = Ci/2K ensure the conservation of probability. Note that, with respect to 
the usual definition, eigenvalues are rescaled by a multiplicative factor —1/2K so that the support of the spectrum 
becomes positive. 



B. Definition of spectral quantities 

We call A e (respectively Wj e ) the eigenvalues (respectively the components of the associated eigenvectors normalized 
to unity) of the Laplacian W s , with e = 1,...,N. Most spectral properties of W s can be obtained through the 
calculation of the resolvent |7| 



G%{X + ie)=[{X + ie)l-W s \ =]T 



N 

Wj,eWk,e 

A — A e + ie 



The mean density of eigenvalues indeed reads 

P(A)= -^e lim + ^ ImG ^ (A + Ze) ' (2) 

where the overbar denotes the average over disorder, that is the random graph Q. 

Another quantity of interest is the autocorrelation function of eigenvectors. The power spectrum |u) e ((9)| 2 of eigen- 
mode e, i.e. the squared modulus of its Fourier transform reads 

N 

\w e {0)\ 2 - Yl w jie w kie exp(2ine(j-k)) . (3) 
j,k=i 

We then define \w(0, A)| 2 as the sum of |u} e (#)| 2 over all u^e lying in the range A < A e < A + dX, divided by the 
number Np(X)dX of such eigenvectors. Once ensemble average is carried out, this power spectrum is simply related 
to the off-diagonal resolvent (|l|) through 



j} E G U X + l£ ) ^M^eu k)) = y°° P ( M ) l ^ 0,t fl , (4) 



where we have assumed that the density of states p is self-averaging |30(] , as confirmed by numerics when the size of 
the sample N becomes large. 

Localization properties of eigenmodes can be studied through the knowledge of inverse participation ratios [^,0 . 
The inverse participation ratios of eigenmodes e, defined as ro* = J^. |w4^| 4 are then averaged over a small energy 
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FIG. 1. The K — 3 ring, a : Relation of dispersion; the plane wave phase <j>(X) is plotted as a function of the eigenvalue 
A(full curve). Non zero stationary eigenvalues are indicated by dashed lines (from left to right : A^, A^ and A^, see text), b : 
Density of states f> TC (A) as a function of A. Van Hove singularities take place at Af, X^, X$ and A^ (from left to right) 

region around A to give u> 4 (A) following the procedure described above. Though inverse participation ratios may 
fluctuate from state to state at a given eigenvalue A, a non vanishing value of the ensemble averaged w 4 (X) gives a 
good indication of the emergence of localization 0] . 

The spectral properties of S will of course reflect its mixed structure. Before turning to numerics, it is therefore 
useful to briefly recall known results on the spectral properties of both TZ and Q separately. 



C. Diffusion on the unidimensional ring 

On the ring, the Laplacian operator W 1 ^ is diagonalized by plane waves whose components at site j read 

( w 4>)j = -j= ex P { 2 ™<f>j) , (5) 
V N 



where the phase <fi assumes multiple values of 1 /N, from up to l — l/N. In the large N limit, <fi becomes a continuous 
variable ranging from < <f> < 1 and the corresponding eigenvalue is given by 



1 - 

A k (0) = 1--]Tcos(2vt^) . (6) 



f=i 



Due to the reflection symmetry <j) — > 1 — 0, we can restrict to <f> belonging to the interval [0; 1/2]. The dispersion 
relation (^) defines the inverse multivalued function (j) n (X) plotted fig la for K — 3. The eigenvalue A is stationary for 
four different wave numbers — 0, 0^ ~ 0.206, cj)^ — 0.354 and =1/2 corresponding to A^ = 0, Xf' ~ 1.439, 
A^ ~ 0.981 and = 4/3. The density of states p n (X) is given by 

and is shown fig lb for K = 3. Van Hove singularities are located at Af, A^, A^ and A^. 

The relation of dispersion may be found back by looking at the autocorrelation function of eigenvectors, or alterna- 
tively at their power spectrum \w(8, A)| 2 , see Section II. B. According to fig la, \w(0, A)| 2 will exhibit one (respectively 
two or three) Dirac peaks over the range < < \ when < A < Xf (respectively Xf < X < Xf or Xf < X < Xf). 
The locations of the peaks coincide with the multivalued function (f R {X) 1 giving back the relation of dispersion (||). The 
study of the power spectrum will turn out to be useful for establishing pseudo-dispersion relations on the small-world 
lattice in presence of disorder. 



D. Diffusion on the random graph 



Diffusion on the random graph has been the object of several analytical and/or numerical studies |jis^4 -28|. Briefly 
speaking, the two main physical features are : 

• The spectral density has a bell-like shape with small oscillating tails (see fig 1 in [p6[ ). The central bell- like part 
corresponds to extended states. 

• The lateral oscillations are accompanied by increase of the inverse participation ratio, which is an indication of 
the localization of eigenmodes near mobility edges. Localized states are centered on geometrical defects, that is 
sites with abnormally large or low connectivities with respect to the average coordination number. 

More quantitatively, defects will pin localized states if their connectivity c deviate from the average coordination 
number q by more than iy/q. Due to the defect mechanism, the localized peak associated to connectivity c has a 
total weight (i.e. integrated density of states) well approximated by the Poisson law p c e~ p /c\. 



E. Perturbation theory and heuristic arguments 

We start by decomposing W s into W n + W Q . One is tempted to treat W s as a small perturbation with respect 
to the ring Laplacian W n when ijCl. The resolvents G and G n for the small- world lattice and the ring obey the 
Dyson equation 

N 

G%(\ + ie) = Gf k (A + ie) + £ Gf e (X + ie) Wf m G s mk {\ + ie) . (8) 

l,m=l 

Using the above equation, the Fourier transform of the self-energy E at energy A and momentum <f> can be computed 
order by order within perturbation theory R. 
At first order, the self-energy reads using 

N 

Ex (A + ie, <t>)=Y, W % Mi = ± , (9) 

3,k=l 

for non zero phases </>. In other terms, the whole spectrum on fig lb would be shifted to the right by q/2K. 
Second order corrections to the self-energy may be written as 

E 2 (A + z £ » = J2 E W % W L Gf t {\ + ie) K)j W m = ^L dn P [ ^ . , (10) 

where the subscript c indicates the connected two-points correlation function of W G . The emergence of an imaginary 
part in the self-energy shows that the resolvent G^ k at energy level A will decrease exponentially with the distance 
\j — k\ over a typical length 

L E (X)= 2 -^q- 1 (p n (X))- 1 . (11) 

The latter may be interpreted as the elastic mean path [§,[3l| travelled by a wave at energy level A between two 
successive scattering events. In other words, the peaks appearing in the power spectrum of the eigenvectors (|^) 
will acquire some finite width ~ 1/Le, see Section IV. D. Le{X) indeed diverges for weak disorder q — > 0. Clearly, 
perturbation theory is sensible as long as the wave length 1/0 TC (A) keeps much smaller than the elastic mean path 
(A), the so-called Ioffe-Regel criterion gj. This condition obviously breaks down at van Hove singularities, see 
Section II. C. Furthermore, at moderately large eigenvalues (A > A^ in the K = 3 case), the existence of several 
distinct branches in the dispersion relation (fig la) makes the density of states higher than at smaller eigenvalues 
(A < Aj") as seen fig lb. We therefore expect from expression (|ll|) that disorder will affect the density of states at 
large eigenvalues more strongly than at small energy levels. 

This reasoning should hold except at vanishing eigenvalues A. One-dimensional model indeed present the peculiarity 
that the density of states diverges at zero energy. More precisely, Le(X) ~ VA as A — > and therefore perturbation 
theory should not be trust below A = O(q). 
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FIG. 2. Results for connectivity parameters K = 3, q = 1. a : density of states p obtained from numerics (full curve) and 
EMA (dotted line) approximation, b : coordination of eigenstate-center c (full curve) and inverse participation ratio w 4 for 
N = 256 (dotted curve), N = 512 (dashed curve) and N = 1000 (dot-dashed curve) averaged over 1000 samples. 

What happens at smaller energy levels? Some heuristic arguments were developed in [T^ ] to answer this equation 
in the random graph case. First of all, calculations of the spectrum on a Cayley tree with fixed connectivity c show 
a gap between the null eigenvalue and the left side of the spectrum support which disappear in the limit c — > 2 
only. Indeed, long enough chain-like structure are capable of exhibiting small eigenvalues of the order of their inverse 
squared length fH| . It is appealing but speculative to think that the small eigenvalues of W s come from rare portions 
of the 1-D ring unaltered by disorder fl^j . This Lifshitz-like argument may be quantified as follows. A chain of M 
sites, giving rise to eigenvalues of the order of 1/M 2 will not be corrupted by the random graph component Q with 
an exponentially small probability scaling as exp(— qM). Summing these contributions to the density of states over 
large M values, we obtain 

p(A)~p*(A)«p(-^) , (12) 

for small A's. Expression ( |l2|) has a nice interpretation in terms of the disordered-induced length = 1/q. The 
eigenmodes of the ring are exponentially attenuated with the ratio of their wave-length over the typical distance 
between two successive random links incoming onto the one-dimensional structure. When A — > 0, the exponential 
term in (|l^) vanishes as exp(— Cq/VX) (where C is a constant) and thus screens the algebric divergence of p n . This 
argument predicts a cross-over between the pure case and the Lifshitz tail at A c . . = 0(q 2 ) corresponding to a density 
Pc.o. = 0(q^ 1 ). We shall see in next Sections that such a cross-over is indeed present. 

III. NUMERICAL DIAGONALIZATION 

In this Section we restrict to K = 3. The corresponding ring has a non trivial relation of dispersion giving rise to 
richer features than for smaller values of the connectivity, e.g. K = 1. In addition, diagonalizations of small- world 
lattices with larger K suffer from bigger finite-size effects and thus reliable results demand a prohibitive computing 
cost. 



A. Description of the eigenvalues spectrum 

We have performed exact numerical diagonalizations of the random Laplacian W s in the case K = 3 for different 
sizes, up to N = 1000 and for 1000 samples. Fig 2a and fig 2b respectively display p(X) and u> 4 (A) for connectivity 
parameters K = 3, q — 1. Main remarks are in order : 
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FIG. 3. Results for connectivity parameters K = 3, q = 5. a : density of states p obtained from numerics (full curve), EMA 
(dotted line) and SDA (dashed line) approximation, b : inverse participation ratio w 4 (full curve) and center coordination c 
for TV = 256 (dotted curve) and N = 512 (dashed curve) averaged over 1000 samples. 



• the density of states at small eigenvalues seems to vanish for A < A_ ~ 0.045 ± 0.005. There is no true gap 
strictly speaking but rather a very small density of states (discussed in Section II. E) that cannot be accounted 
for by numerics due to the statistical shortfall of eigenvalues. In the following, the [0; A_] range will be referred 
to as a "pseudo-gap" . 

• the density of states exhibit a maximum on the left side of the spectrum at A c o . ~ 0.11, slightly above A_. 

• the central part of the spectrum (A c . . < A < A+) has a smooth shape and corresponds to extended states. The 
general form of p(X) is reminiscent of the density of states for the ring fig lb but divergences have been smeared 
out by disorder M. Notice that in the range A c . . < A < 1, the density of states is in quantitative agreement 
with the ring spectrum, compare fig lb and fig 2a. For increasing sizes N and at fixed A, w A vanishes as 1/N 
and the breakdown of this scaling identifies the upper mobility edge : A+ ~ 2.1 ± 0.1. 

• on the right side of the spectrum, the density exhibits successive regular peaks and the eigenstates become 
localized, see fig 2b. 

Similar results are shown fig 3a-b for K = 3 and a larger random graph component q = 5. In this case, we have 
found A_ ~ 0.4 ± 0.05, A c . . ~ 0.7 ± 0.05 and A + ~ 3.3 ± 0.1. As expected, localization effects are stronger and the 
pseudo-gap becomes wider. 



B. Relations of dispersion for extended states 



Following Section II. B, we have computed the power spectra \w(9, A)| 2 for different values of A and connectivity 
parameters K = 3, q = 1 and K = 3, q = 5, see fig 4 and fig 5 respectively. Depending on the eigenvalue A, the power 
spectrum is mainly composed of different peaks whose number varies from two to six as in the ring case, see fig la. 
However, these peaks have now non-zero widths, equal to the inverse autocorrelation lengths of the eigenvectors. A 
good fit on the range < 9 < 1 consists in a superposition of 2 m Lorentzian spectra, see fig 4 and fig 5, 

m 

Peaks are labelled so that < <j)\ < ... < <j> m . Close to the pseudo-gap edge A_, m equals two and the weight 
W2 of the secondary peak gets much smaller than wi 7 see fig 4a and fig 5a. As a result, the estimates for §i,<Ji 



W( 



Att 2 {9- fa) 2 + aj 4tt 2 (1- 9 - fa) 2 + , 



(13) 



G 




FIG. 4. Power spectrum \w(0, A)| 2 of the eigenvectors as a function of the phase 9 for K — 3, q = 1 at different energy levels. 
The dotted lines correspond to the Lorentzian fit ([13[) of Section III.B. a : A = 0.8. b : A = 1.2. c : A = 1.8. 




FIG. 5. Same as fig 4 but with connectivity parameters K — 3, q = 5. a: A = 0.9. b : A = 1.7. c : A = 2.5. 



becomes less accurate. Moreover, for large eigenvalues, the Lorentzian form (13) does not fit numerical data two well. 
Consequently, the error bars on fa, og increase when A reaches A + . We have plotted fig 6 and fig 7 the wave numbers 
fa and the widths at obtained from the numerical data vs. the eigenvalue A for q = 1 and q — 5 respectively. 

At first sight, the wave numbers <pg are much less affected by disorder than the overall features of the density of 
states, e.g the edges A_ and A + and of course the inverse lengths ag (that are equal to zero in the ring case). Indeed, 
the stationary points in the wave numbers curves fig 6a and fig 7a, that is the merging points between fa, fa and 
4>2, fa precisely coincide with the ring values fa* 1 and fa% of Section II. C. The main qualitative difference between the 
above figures and fig la lies in the appearance of a narrow area surrounded by the fa and fa lines, centered around 
<j>2 and extending down to A = A_. Thus, contrary to the ring case, the power spectrum of these eigenstates displays 
more than one peak (even if the heights of the secondary peaks are small as mentioned above) . 

When A lies in the range A c . . < A < 1, o\ gets very small, see fig 6b and fig 7b. In presence of the disorder, these 
eigenmodes have conserved an overall quasi-plane wave form. Such a behaviour is indeed expected from Section II. E. 
This result explains the quantitative agreement between the density of states of the ring and the numerical spectrum 
for the small- world lattice observed in Section III. A in the intermediate range A c . . < A < 1. 

To end with, let us stress that formula (13) allows us to define some pseudo relations of dispersion characterizing 
the gross form of eigenmodes at energy A through some wave numbers fa and inverse autocorrelation lengths og. 



C. Localized eigenstates 



Following Section II. D, we have measured for each eigenvector wi the connectivity eg of its center, that is the site io 
with maximum component |wj 0) e| |26| - The mean connectivity c(A) of the centers of eigenvectors having eigenvalue A 
is plotted fig. 2b. It is a smooth monotonous function of A in the central part of the spectrum. In the localized region, 
c(A) is constant over a given peak and integer- valued (for c > c+ = 10 on the right side); the center connectivity 
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FIG. 7. Same as fig 6 but with connectivity parameters K = 3, q = 5. 



8 



c 


1 1 
11 


1 O 


1 Q 


1 A 
14 


^NUM 


2.175 


2.315 


2.460 


2.615 


X SDA 


2.195 


2.325 


2.469 


2.619 


p NUM 


0.60 


0.78 


0.90 


0.44 


pBDA 


0.550 


0.686 


0.755 


0.798 



TABLE I. Case K = 3, q = 1 : Weights p (i.e. integrated density of eigenvalues belonging) of the localized peaks (divided 
by the Poisson factor e~ q q c /c'\ with c' = c — 2K) and corresponding eigenvalues A for defect connectivities c. Results are 
obtained from numerical simulations (NUM) and SDA. The accuracy of numerical values is discussed in Section III.C. 



c 


18 


19 


20 


21 


22 


^JV U M 


3.425 


3.545 


3.685 


3.850 


3.990 


X SDA 


3.433 


3.556 


3.696 


3.841 


3.990 


p NUM 


0.60 


0.62 


0.79 


0.77 


0.60 


p SDA 


0.574 


0.696 


0.762 


0.802 


0.831 



TABLE II. Same as Table I but with connectivity parameters K — 3, q = 5. 



abruptly jumps when A crosses the borders between peaks. 

Table I and II list the weights p^ UM of the localized peaks related to connectivities c, that is the integrated densities 
of eigenvalues belonging to them. The corresponding eigenvalues A^ are also given. They are measured at the top 
of the peaks with absolute error ±0.005 whereas the relative error on p^ UM is about 20%. The reliability of numerical 
results suffers from two effects. First, as A decreases, localized peaks get closer and closer to the flank of the extended 
spectrum and the determination of the top of the peaks \^ UM becomes less and less accurate. Secondly, for large 
eigenvalues, the corresponding statistical events are rare and fluctuate drastically from sample to sample leading to 
a poor accuracy on p^ UM ■ Note however that p^ UM is of the same order of magnitude (but smaller as it should be 
due to the screening effect exposed in [|(|) as the fraction of sites having c = 2K + c' neighbors, with c' given by a 
Poisson law of parameter q. 

Contrary to the random graph case J2(|, localized peaks are not found on the left side of the spectrum. Indeed, the 
minimal connectivity of a site in the small-world case is c = IK and defects with low connectivities with respect to 
the average coordination IK + q can be seen at very strong disorder only. 



IV. THEORETICAL FRAMEWORK 



A. The average over the disorder 



Once averaged over the disorded, the resolvent (|l|) can be written as the propagator of a replicated Gaussian field 
theory § 



Gf k (X + ie) = lirn -- J dxi(xj.a!k) exp ( -(A + 



^£3 + l£w«(*-2<) 2 



/ J I 



(14) 



Replicated fields Xi are n-dimensional vector fields attached to each site i. 

The Laplacian on the ring 1Z is non fluctuating and we have to perform the average over the random graph Q only. 
Due to the building process described in Section II. A, the averages over different pairs of points are decoupled and 
we find 



exp r-^wf^-xA =n 

V i<£ / Ki 

~ exp 



1 exp \Xj — xp 



qN a 
2 N 



^exp 



AK 



(xi - x e y 



(15) 



to largest order in N. An important property of the last term in (|l_5j) is that all pair of sites i, i interact together with 
the same coupling. Therefore, introducing the density p(x) of replicated vectors Xi equal to x, 
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N n 
i—l a—1 

by means of a functional Lagrange multiplier u?(x), we obtain 



(16) 



exp 



|E^ 



a - 



i<£ 



= / V ' p{x)T>u>{x) exp 



iN 
2~ 

qN_ 
2 



dx u)(x)p(x) + - uj(xi) — 



N 



qN 
~2 



dxdy p(x)p(y) exp ( ~-^{x - yf 



(17) 



where the above expression includes a double path-integral over p and u> functions. We shall come back in Section 
IV. C to the significance of the latters. 

We can now take into account the contributions coming from the diagonal term and from W n in (]l4|). The n th 
moment of the partition function Z of the replicated Gaussian field theory introduced in ( |l4| ) 



may be written as 



Z n = I Vp(x)Vuj(x) expi-NTip,^)) 



with 



dx uj{x)p{x) 



dxdy p(x)p(y) 



exp 



(x-y) : 



K 



log A m {oj) 



(18) 



(19) 



(20) 



In the above expression, Aj\/(w) denotes the largest eigenvalue of a transfert matrix T we shall define in next section 
and which depends on the function lj. Note that (^o|) is exact if N is a multiple of K. 

The analogy with the theory of polymers initiated by Edwards |[6]|l7]] is clear once the p variables have been 
integrated out in (|l9|). We are then left with a one-dimensional model of free-energy density — log Aju(w)/if in 
presence of a random potential u distributed according to a Gaussian law. For the small-world lattice, this is no 
approximation since the virial expansion includes a unique corrective term to the non-interacting case, see Section 
IV.C. 



B. Transfer matrix 

Each element of the transfer matrix T arising in the computation of the free-energy functional ( po| ) is labelled by 
IK replicated vectors x\, a?2, • • • , %k an d VuV2, • ■ ■ ,Vk- From ( p^ ) and (p"5|), T reads 

!K . K 

-(A + it) ^[z/ + yf\ + - ^2[v(xj) + uj(yj)] 
3=1 3=1 

-4k E [ft-^ + t^^-^E&Vffi) 2 • (2i) 

i<j<e<K j=i i=\ J 

Due to the last term in (|2l]), T is not a symmetric matrix. However, defining its "mirror" matrix T# through 

T*[xi,x 2 , ■ ■ .,x K ;yi,y 2 , ...,y K ]= T[y K) y K -i, . . . ,yi;xi,x 2 , ...,%] , (22) 

it appears from (|2l]) that T = Let us now define the mirror v# of a vector v through 

v*[yi,m, ■ ■ ■ ,Vk] = v[y K ,yK-i, ■ ■ ■ ,yi] ■ (23) 

Then, it can be easily checked that the left eigenvectors of T are the mirrors of its right eigenvectors. Consequently, 
in the large limit, the matrix entries of T N simplify to 

(T N )[xi,x 2 , . . .,x K ;m,y2, . . . ~ (Am)* ^[2/1,2/2, ■ ■ ■ ,yK)vM[xi,x 2 , ■ ■ ■ ,xk] , (24) 
where vm is the maximal eigenvector of T associated to Am- 
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C. Meaning of the order parameters p and u 



In the large N limit, the n th moment of the partition function ( )lS| , |l9| ) can be calculated using the saddle-point 
method. The optimization of the free-energy functional J- (^) gives the values of p(x) and u>(x), as well as 
vm[xi, x 2 , • ■ ■ , xk] at saddle-point. 

The significance of p is obvious from definition fll6| ); this is the probability distribution of replicated vectors Xj. 
Therefore, p gives access to the average resolvent (|l4|) 



lim 

n— >o n 



dx p(x) x 



(25) 



and to the density of states through (||) . 

More generally, within the transfer matrix formalism developed in Section IV. B, the joint probability distribution 
of replicated vectors x*j,xV|-ij . . . ,x*i+K-i is related to the maximal eigenvector vm of T. Normalizing vm to unity, 
this joint probability reads (B4|) 



Pjoint(Xi,X i+1 , . . . ,X l+K -l) = V M [Xi,X i+1 , . . . ,S i+ K-l] vff[Xi,Xi+t, . . . ,X i+ K-l] ■ (26) 

In particular, the one point density p is found back when all but one vectors are integrated out in the above formula 



p(x) = J dx 2 ... dx K v M [x, x 2 ,..., x K ] v%[x, x 2 , ■ ■ ■ , x K ] 



(27) 



Identity (^7|) is precisely the stationary condition of T (20) with respect to uj. 

The meaning of ui is less straightforward but emerges from the extremization condition of T with respect to p, 



q / dy p(y) 



cx p(-^( f -y) 2 ) - 1 



(28) 



The above equation is stronglyreminiscent of the first correction to the logarithm of the density obtained from Mayer's 
expansion of interacting gas |3^] . Indeed, this close relationship has been exposed in |2(| (see also [[33| for the similar 
case of finite-connectivity spin-glasses) . The replicated field Xi attached to site i may be interpreted as the position of 
particle i in a n— dimensional abstract space. Particles see an external harmonic potential — |(A + ie)x 2 and interact 
with each other through a two-body quadratic energy — f (x} — Xj) 2 , see (|l8|). In the absence of any W s matrix, 
the density of a non-interacting gas is recovered 



p (x) = exp ( ^(A + ie)x 2 



(29) 



obtained from equation (|l8|). The first virial coefficient, that is the first correction to \og p(x) in the expansion in 
powers of p (and not in powers of the fugacity po) is precisely given by the r.h.s of ( |28| ) p6] , p2| . It turns out that 
Mayer's expansion for random graph do not include higher order terms, as is expected from the similarity between 
random graph and trees. 

Therefore, —^lo(x) may be understood as the effective potential due to interactions with other particles seen by a 
particle located in x. This interpretation will be useful in deriving the single defect approximation of Section VI. 



D. Autocorrelation function of the eigenvectors 

We have seen in Section II. B that the autocorrelation functions of eigenvectors are related to the resolvent G^ fe , 
see identity (Q). By definition, the latter is the average scalar product of replicated vectors Xj and Xk ( |l4] ) with the 
Gaussian measure (|l8|). Within the transfer matrix formalism we have introduced in Sections IV. B and IV. C, this 
mean dot product can be computed from the knowledge of all eigenvectors vt and eigenvalues Ae of T. We decompose 
the distance d = \j — k\ between sites j and k as d — do + K d\ with < do < K — 1 and d\ integer- valued. Then the 
average resolvent reads 

°° / A \ di 

G%(\ + te) = to*-i y E(v*.Xlv t ) (vf.X± +do .v M ) ( J L ) • (30) 
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In the above equation, the maximal eigenvalue of T coincides with £ = 1 (Ai = Am) and increasingly excited states 
correspond to I > 2. The X£ operator measures the value of the a th component of the replicated field x appearing at 
the b th place in ve (with 1 < b < K — 1) and (vfj.X.Vi) denotes the matrix element of between states vm and vg, 

{v%.X%.Vi) = J dxi...dx K x% vfj[xi,X2,...,x K ]vi[xi,X2,...,x K } ■ (31) 

For small distances d< K — 1, d\=Q and the sum over eigenstates in formula (|3(]) can be carried out to obtain 

GfJX+Jej = lim -i {v*.X\ X\ +d .v M ), Q<d<K-l. (32) 

Using identity the above equation gives back the d — resolvent (|25|). 

Expression (|30|) may also be used to compute the autocorrelation function of eigenvectors at large distance d, that 
is when the sum in ( J3 0|) is dominated by the 1 = 2 contribution. The ratio r = A2 / Am whose modulus is by definition 
smaller than (or equal to) unity can be written as 

r = cxp {-Ko + i2irK<j)) , (33) 

with a > 0. Therefore, at large distances d, the average resolvent scales as 

Gf k (X + ie) ~ exp(-cr d + ilv^ d) , (34) 

where both a and cj> depend on the energy level A. Consequently, a may be seen as the inverse autocorrelation length 
of the eigenmodes associated to level A and 4> to the typical wave number of the latters. We shall explicitely compute 
the inverse length a and wave number <f> in Section V.D within the effective medium approximation and compare 
them to the numerical results of Section III.B. 



V. EFFECTIVE MEDIUM APPROXIMATION 



The exact extrcmization of free-energy function (|20|) is an awkward task. We start by solving this problem using 
an effective medium approximation (EMA) which gives a good description of the extended part of the spectrum. 



A. Presentation of the approximation 

The starting point of EMA is the exact expression of the density p obtained from (|18| 



^) = lgC,exp( 2[(A+ .^_^ ] _ 1 ) , (35) 

where the Cj are normalization constants going to unity as n vanishes. 

In the extended part of the spectrum, we expect all matrix elements appearing in ( pq ) to be of the same order of 
magnitude and thus p(x) to be roughly Gaussian. We therefore assume the following Ansatz for the density 



^^)=(2^i 9 EMA )--exp\^ MJ ) . (36) 

In addition, the effective potential u> we choose to be harmonic, see Section IV. C, 

u EMA (x)=g EMA x 2 . (37) 



i x 



2 



EMA is therefore implemented by inserting the Gaussian Ansatz (|36|,|37|) into functional T. Due to the choice (|3 
for lo, the transfer matrix T simplifies to the tensorial product of n identical matrices 



K 

T EMA [ Xl ,x 2 , x K ; yi ,y 2 , ...,y K ]= exp { -{g EMA + A + ie) + y 



3=1 

i<j<e<K j=i i=\ 
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The resulting free-energy per replica reads in the limit n — > 0, 

1 a ( n EMA \ 1 



2" 3 4°V K J K 

where Af { MA {g EMA ) is the largest eigenvalue of T EMA @. 

B. Self-consistent equations 

To compute A|J f M ' 4 (g- EM " 4 ), we take advantage of the Gaussian structure of T EMA and look for a maximal eigen- 
vector of the form 



v M [xi,x 2 ,...,x K \ =v .exp 



O X! V " r X Q Xr ) ' ( 40 ) 

\ q,r=l J 



where vq is a normalization constant. The self-consistent equations fulfilled by the K x K matrix V are immediately 
obtained and reads 

V ' 8=1 t=l 

CV = V + 5 V (ff £MA + A + — - i) + ± , (41) 



where S denotes the Kronecker symbol. The saddle-point equations of F EMA with respect to g EMA and g EMA 
respectively read 

a EMA - 1 (42) 

9 2(g EMA - K) [ ' 



EM A = 1 T 

y K 



V + V 



(43) 



The saddle-point equation (27) gives back equation (|43|) while identity ( |42| ) differs from (|28|) since the EMA order 
parameters (36 3^) are not exact solutions of the extremization conditions over T . 

We have solved numerically the above self-consistency equations using the followingprocedure. For given K ', q and 
A, we start with an arbitrary value of g EMA , e.g. null, and compute g EMA from ([42]) . We then seek for matrices 
V and U fulfilling ([i"l|). Among all possible solutions, we select the matrix V having all eigenvalues with positive 
imaginary parts so that v M MA (^5|) is bounded for large arguments x. The new value of g EMA is computed from jiij). 
The whole process is iterated until a fixed value of g EMA has been found. The density of states is then equal to the 
imaginary part of g EMA divided by — n, see in order (Bff), (p5[) and (0). 



C. Spectrum, mobility edge and pseudo-gap 

The density of states obtained for K = 3, q — 1 and q — 5 are shown fig 2a and fig 3a respectively. The overall 
form of the EMA spectrum coincides with the numerics. The quantitative agreement is of course less precise at large 
disorder (q = 5). We shall see in Section VI. C how the theoretical spectrum may be improved using SDA. 

On the right side of the spectrum, the mobility edge is accurately estimated : X EMA ~ 2.13 for q = 1 and 
\EMA ^ 3 f or q — 5 ; compare to Section III. A. Nevertheless, as expected from the presentation of EMA in Section 
V.A, the oscillations of p{\) corresponding to non-extended states are not captured by the approximation. 

A similar situation takes place at small energy levels. EMA predicts a gap whose order of magnitude (\^_ MA ~ 0.064 
for q = 1 and \ EMA ~ 0.46 for q = 5) coincides with the findings of Section III. A. However, from a qualitative point 
of view, EMA wr ongly p redicts a vanishing density of state below \ EMA . This artifact is well known in the case of 
the Bethe lattice |T^,|l2],g9[ for which EMA becomes exact. 

We have analytically checked the scaling behaviour discussed in Section II. E in the simplest case K = 1. In the 
small q, A region the EMA resolvent is obtained from Section V.B, 
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The gap width is therefore given by \^_ MA ~ q 2 /8, in agreement with the heuristic arguments of Section II. E 
indicating that the size of the pseudo-gap scales as 0(q 2 ) for weak disorder. Furthermore, the cross-over between 
"one-dimensional" behaviour and rare events take place when the imaginary part of the resolvent ( pl4| ) reaches its 
maximum, that is at A c . . — 9 2 /4. The corresponding density of states scales as p c . . ~ 1/q, in quantitative agreement 
with the findings of Section II. E. Notice also that the value of the ratio A c . ./A_ ~ 2 is supported by the numerical 
results for K = 3, see Section III. A. 



D. Relations of dispersion 

We have shown in Section IV. D that the autocorrelation function of eigenvectors is accessible from the second largest 
eigenvalue A2 of T (^l|). Due to the similarity of T EMA ([38]) with the evolution operator of a quantum oscillator, we 
look for an eigenvector of the following form 



,EMA 



(l K \ K 

lrX r , (45) 

\ q,r— 1 / r— 1 



where V\ is a normalization constant and V qr is the quadratic form appearing in the maximal eigenvector vm and 
self-consistently defined through ( pll"| ) . An elementary calculation shows that the above Ansatz indeed corresponds to 
an eigenvector of T EMA if and only if the K coefficients 7 r arising in (|45|) fulfill the eigensystem 



K K 

^r = -^J2J2( U ~^^ > ( 46 ) 
q—r s—1 

where U has been defined in © and r = A 2 /A M . We are thus left with the diagonalization of a K x K matrix 
whose complex eigenvalues 73, t = 1,...,K can be written as (Q), 

T( = exp (-Kae + i2-n-K4>i) . (47) 

In doing so, we have at our disposal the first K excited states of T and not only the first one. This is a peculiarity 
of the EMA scheme, that was not expected from Section IV. D. The average resolvent therefore obeys the asymptotic 
scaling ( ^i| ) 

K 



Gf h (\ + ie)~Yl Ce ex P(~^ d + i27T( t>e d ) ■ ( 48 ) 
1=1 

This large distance behaviour strengthens the data modelling (|l^) used in Section III.C. Indeed, if relation 
held for all d's, the power spectrum of the eigenvectors \w{0, A)| 2 would be well approximated by a superposition of 
Lorentzian factors from equation (Q). 

The wave numbers <f>i as well as the inverse autocorrelation lengths ui are plotted fig 6 and fig 7 vs. the eigenvalue 
A for K = 3, q = 1 and K = 3, q = 5 respectively. The agreement between the EMA prediction and the numerical 
results of Section III.B is good. Notice that numerical points for the wave numbers (fix and 4>2 merge and lie between 
the theoretical branches when A approaches A+. This phenomenon is simply due to the coalescence of the first two 
peaks in ( p"3] ) when \<p2 — <j>\\ — <7i/27r arising at such energy levels. 

In the extended part of the spectrum, EMA can therefore be used to obtain some accurate effective relations of 
dispersion for the eigcnmodcs of the Laplacian operator, even at strong disorder. 



VI. SINGLE DEFECT APPROXIMATION 



A. Principle of the approximation 

The physical grounds of the single defect approximation have been exposed in [ p6[ . The basic idea is to treat 
exactly, via the stationarity equations ( p7| , p8| ) the interactions of a single site, i.e. the defect with its surrounding 
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neig hbors belonging to the effective medium |29| taking into account the fluctuations of the number of these neighbors 
|26| . This approximation scheme amounts to perform one iteration of the extremization condition on the free-energy 
functional T from the EMA solution. Let us see how SDA works for the small-world problem. 

First, we consider the saddle-point equation for the effective potential d2§|). When inserting the EMA density (|3£ 
in the r.h.s. of (Eq), we obtain the SDA effective potential 



— ui 
2 



SDA 



(x) = q 



exp 



EMA 



2K) 



1 



(49) 



up to irrelevant terms when n 
To calculate the SDA density p SDA 



0. This expression differs from the best quadratic potential obtained from EMA 



the stationarity identity (^7|) requires the eigenvector of the transfer matrix (|2 



in presence of the SDA potential (49). Such a calculation cannot be done analytically. Remembering the interpretation 
of u> exposed in Section IV. D, we may however estimate the ratio of the probabilities that a particle be located in the 
n-dimensional abstract space at position x within each scheme of approximation by 

W ~, „„„ , .SDA 



„EMA 



(X) 



exp 



2 W 



exp 



EMA 



2K 



-EMA 



(50) 



As in the simpler case of the random graph |p6| , the anharmonicity of the effective potential uj sda ( fl9|) reflects the 
local fluctuations of connectivity. Using the above expression for p SDA , we can now study the localized states induced 
by the geometrical defects we have observed in Section III.C. 



B. Localized eigenvalues 



The second moment of the field x in the SDA scheme reads, see (E 



..SDA 



n 
e~9 



= lim dx p bUA (x) x 

oo 

£ 

c=0 



-I 



n EMA 



n EMA 



— 2K 



)E ma\ 



(51) 



Above the mobility edge X^ MA , the EMA density of states vanishes. Thus, g EMA (X) and g bUA (X) are real numbers. 
However, the poles of g SDA located on the real axis can give rise to Dirac peaks in the density of states p6fl . Consider 
indeed such a singularity taking place at X S ing > \+ MA . Close to the singularity, we may write 



..SDA 



(A) 



A- A 



(52) 



sing 



and obtain a contribution to the density of states equal to p S mg 8(X — X S i ng ) using (|2j). As seen from expression (|5^ 
the poles of g SDA form a discrete sequence of energy levels Af DA labelled by the integer number c. The values of 
Af DA and of the corresponding weights p^ DA arc listed Table 1 and Table 2 for K = 3, q = 1 and K — 3, q = 5 
respectively. The agreement with the numerical values of Section III. A is good (see that Section for a discussion of the 
reliability of numerical measures). Note that c coincides with the coordination number of the centers of eigenmodes 
measured in Section III.C. This strongly supports the statement that defects are responsible for localization |2q|. 



C. Extended spectrum at strong disorder 

The EMA density of states for K — 3, q — 5 displayed fig 3a shows some discrepancies with the numerics from the 
maximum density located at A ~ 1.5 — 2.0 up to the mobility edge. We have computed the SDA density of states in 



the extended region from equations (51) and (|2|). The result is shown fig 3a and agrees in a much better way with 
numerics. 

For completeness, we have also compute the SDA density of states at low disorder K = 3, q = 1 in the extended 
region. The resulting spectrum can hardly be distinguished from the EMA prediction, and is therefore in very good 
agreement with numerics. 
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VII. SUMMARY AND CONCLUSION 



In this work, the spectral properties of the Laplacian operator on mixed lattices, composed of a random graph 
structure superimposed to a unidimensional ring have been investigated. Using some approximation schemes, we have 
been able to reproduce the numerical density of states with a good agreement and to obtain a quantitative description 
of the localized eigenstates. From this point of view, we have shown that the single defect approximation introduced 
in the simpler context of the random graph was also capable of giving reliable and precise estimates of the localized 
energies in presence of a finite-dimensional underlying geometry. 

One of the most interesting aspects of the present work is the existence of pseudo dispersion relations for disordered 
and thus non periodic lattices. By studying the autocorrelation functions of eigenmodes, or more precisely their power 
spectra we have shown that some features of the pure lattice remain present even at strong disorder. It is encouraging 
to notice that the effective medium approximation already provides a very good description of the power spectra of 
the eigenmodes. To what extent this result holds in higher dimensional models is of course questionable. 

As emphasized in the introduction, the small-world architecture may be of particular relevance to the field of 
polymers and more generally chain-like systems with three-dimensional interactions, e.g. proteins, DNA, ... With 
respect to the diffusion problem studied here, realistic models of polymers or biological molecules have to take into 
account the internal degrees of freedom of the monomers/constituents and their non-trivial interactions. We hope 
nevertheless that the technical approach developed in the present paper will be of use in these more complex cases. In 
particular, the knowledge of the excited states of the self-consistent transfer matrix should give access to the effective 
dispersion relations of the eigenmodes. The latters are of primary importance to understand the dynamical properties 
of such molecules 34 , |j| , with possible applications to the relaxation of single DNA molecules |3q] . 
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